Boltzmann legacy and wealth distribution 
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We briefly review results on nonlinear kinetic equation of Boltzmann type which describe the evolution of 
wealth in a simple agents market. The mathematical structure of the underlying kinetic equations allows to use 
well-known techniques of wide use in kinetic theory of rarefied gases to obtain information on the process of 
relaxation to a stationary profile, as well as to identify simple interaction rules which are responsible of the 
formation of Pareto tails. 



I. INTRODUCTION 

In kinetic theory of rarefied gas, the evolution of the phase 
space density is described by the Boltzmann equation Illll[l2ll 

^/(r;x,y) = -v.V,/(f;x,y) + e(/(r;x,v)). 

This equation contains terms accounting for the two ways that 
the density can change. The 

-yV.,/(f;x,v) 

term represents the effects of streaming; that is, the motion 

xo 1-^ xo + (t - to)vQ vqh^vo (1.2) 

of molecules between collisions. The Q{f{x, v, t)) term repre- 
sents the effects of binary collisions and describes relaxation 
to the local Maxwellian equilibrium rlll [l2ll . 

The collision term Q accounts for all kinematically possi- 
ble (those that conserve both momentum and energy) binary 
collisions. The post-collisional velocities v* and w* are given 
by 



V* = — (v + w+ Iv- 
2^ 



- w n 



w 



-(v + vv— |v — wjii), 

(1.1) 



in which n is a unit vector 
For Maxwellian molecules 
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B((v^ 
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■n)[f{v*)g{w*)]dwdn-pf{v). 

(1.2) 

B{-) is a measure of the collision frequency, dn denotes the 
normalized surface measure on the unit sphere S^, and 



p= j f{v)dv 



The Boltzmann description relies on a main ansatz- While the 
number of gas molecules is enormously high so that collisions 
between molecules do not retain memory, in view of rarefac- 
tion only binary collisions are important for the evolution. 

Starting from the pioneering works of Mandelbrot |28], it 
is now commonly accepted by the kinetic community that 
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in many aspects a trading market composed of a sufficiently 
large number of agents can be described using the laws of sta- 
tistical mechanics, just like a rarefied gas, composed of many 
interacting particles. In fact, there is an almost literal transla- 
tion of concepts: molecules are identified with the agents, the 
particles' energy correspond to the agents' wealth, and binary 
collisions translate into trade interactions. This modelling is 
clearly rather ad hoc, but if one is willing to accept the pro- 
posed analogies between trading agents and colliding parti- 
cles, then various well established methods from kinetic the- 
ory and statistical physics are ready for application to the field 
of economy. Most notably, the numerous tools originally de- 
vised for the study of the energy distribution in a rarefied gas 
can now be used to analyze wealth distributions. In this way, 
the kinetic description of market models via a Boltzmann- 
type equation provides one possible explanation for the de- 
velopment of universal profiles in wealth distributions of real 
economies. 

The distribution of wealth in a market economy is nowa- 
days of great interest. Among other approaches, the de- 
scription of market models via kinetic equations is a fertile 
ground for research, well documented by numerous contribu- 
tions in the recent books I1 3[[l 4l USl [33l [3411 . or by the intro- 
ductory articles 1 15, 24, 25', [Sg]. The bridge between Maxwell 
molecules and nonconservative economies has been outlined 
by F. Slanina |32ll. In this paper, a clear parallelism between 
the evolution of wealth in a simple economy and the evolu- 
tion of the particle density in a one-dimensional dissipative 
gas has been established. This paper motivated to eventually 
adapt more and more of the ideas, which have been developed 
in the studies of dissipative Maxwell gases, to the economic 
framework. 

It should be emphasized, however, that there are substantial 
differences between the collision mechanism for classical gas 
molecules and the modelling of trade interactions. In the new 
framework, interactions can lack the usual microscopic con- 
servation laws for (the analogues of) momentum and energy; 
moreover, random effects play a crucial role. In fact, the key 
step in establishing a reasonable kinetic market model is the 
definition of sensible rules on the microscopic level, i.e., the 
prescription of how wealth is exchanged in trades. Such rules 
are usually derived from plausible assumptions in an ad hoc 
manner. (This is clearly in contrast to the original Boltzmann 
equation, where the microscopic collisions are governed by 
the laws ( II. 11 1.) 

Nevertheless, strong analogies remain. Maybe the most im- 



portant one is related to the fact that, as it happens in the 
classical Boltzmann equation iflU [l2ll . where relaxation to 
Maxwellian equilibrium is shown to be a universal behavior 
of the solution, here the corresponding output of the model 
are the macroscopic statistics of the wealth distribution in 
the society, to which the solution is shown to relax. Mostly 
important, while relaxation to equilibrium in the Boltzmann 
equation is achieved by looking at the monotonicity of the en- 
tropy functional, relaxation to the steady wealth profile can be 
achieved by looking here at the monotonicity of new convex 
functionals. 

The comparison of this output with realistic data is up to 
now the only means to evaluate — a posteriori — the quality 
of a proposed model. For instance, it is commonly accepted 
that the wealth distribution should approach a stationary (or, 
in general, a self-similar) profile for large times, and that the 
latter should exhibit a Pareto tail. Such overpopulated tails are 
a manifestation of the existence of an upper class of very rich 
agents, i.e. an indication of an unequal distribution of wealth. 
The various articles in lUSIl provide an overview over histori- 
cal and recent studies on the shape of wealth distributions; see 
also 1 15] for a collection of relevant references. 

In general, the richness of the steady states for kinetic mar- 
ket models is another remarkable difference to the theory of 
Maxwell molecules |8]. While the Maxwell distribution is the 
universal steady profile for the velocity distribution of molec- 
ular gases, the stationary profiles for wealth can be manifold, 
and are in general not explicitly known analytically. In fact, 
they depend heavily on the precise form of the microscopic 
modelling of trade interactions. Consequently, in investiga- 
tions of the large-time behavior of the wealth distribution, one 
is typically limited to describe a few analytically accessible 
properties (e.g. moments and smoothness) of the latter 

A variety of models has been proposed and numerically 
studied in view of the relation between parameters in the mi- 
croscopic rules and the resulting macroscopic statistics. The 
features typically incorporated in kinetic trade models are sav- 
ing effects and randomness. Saving means that each agent is 
guaranteed to retain at least a certain minimal fraction of his 
initial wealth at the end of the trade. This concept has prob- 
ably first been introduced in |13] (see also LUiZQJ), where a 
fixed saving rate for all agents has been proposed. Random- 
ness means that the amount of wealth changing hands is non- 
deterministic. Among others, this idea has been developed in 
ifign . in order to include the effects of a risky market. Depend- 
ing upon the specific choice of the saving mechanism and the 
stochastic nature of the trades, the studied systems produce 
wealth curves with the desired Pareto tail — or not. 

In this short review, we describe the Boltzmann-like ap- 
proach to wealth distribution, and compare a selection of re- 
cently developed models. We will mainly treat two different 
types of interactions. The first type is such that the binary 
trade is microscopically conservative, while the second type 
is such that the binary trade is conservative in the statistical 
mean only. In these situations, the mean wealth in the model 
Boltzmann equation is preserved, and one expects the forma- 
tion of a stationary profile. 

In the class of pointwise conservative trades, we focus on 



the model designed by Chakraborti and Chakrabarti fl3j, and 
on variants of it. In the class of conservative in the statis- 
tical mean, we focus on models with risky investments, origi- 
nally introduced by Cordier, Pareschi and the author [19]. The 
applied analytical techniques, however, easily generalize to a 
broader class of conservative economic games. These tech- 
niques have been applied in the current mathematical litera- 
ture ][l9, 21-23, 29-31], where kinetic econophysics has been 
treated in the framework of Maxwell-type molecules. The in- 
terest reader, who wishes to obtain a deeper understanding of 
the mathematical roots (and possible extensions) of the ap- 
plied tools, is referred e.g. to Ill0ll35[l . 

The kinetic approach is complementary to the numerous 
theoretical and numerical studies that can be found in the re- 
cent physics literature on the subject, from which it differs in 
several subtle points. In our point of view, the evolution of 
wealth dendity is entirely based on the spatially homogeneous 
Boltzmann equation associated to the microscopic trade rules 
of the respective model. Agents on the market are treated as 
a continuum, just like molecules in classical gas dynamics. 
Not only does this approach constitute the most natural gen- 
eralization of the classical ideas to econophysics. But more- 
over, it clarifies that certain peculiar observations made in en- 
sembles of finitely many agents and in numerical experiments 
(like the apparent creation of steady distributions of infinite 
average wealth, e.g. lllSMlTll ) are genuine ^«/fe size effects. 

For the sake of completeness, a comment on the justifica- 
tion of kinetic market models is in place. The socio-economic 
behavior of a (real) population of agents is extremely com- 
plex. Apart from elements from mathematics and economics, 
a sound description — if one at all exists — would neces- 
sarily need contributions from various other fields, including 
psychology. Clearly, the available mathematical models are 
too simple to even pretend to reflect the real situation. How- 
ever, the idea to describe economic trades in terms of a kinetic 
equation gives rise to a variety of challenging mathematical 
problems, both from the theoretical and numerical point of 
view. In particular, it is remarkable that this class of simple 
models possesses such a wide spectrum of possible equilib- 
ria (some of which indeed resemble realistic wealth distribu- 
tions). Moreover, kinetic market models are extremely flex- 
ible with respect to the introduction of additional effects. In 
this way, the described models should be considered as ba- 
sic building blocks, that can easily be combined, adapted and 
improved. 



II. BOLTZMANN MODELS FOR WEALTH 

We consider a class of models in which agents are indis- 
tinguishable. Then, an agent's "state" at any instant of time 
f > is completely characterized by his current wealth w > 0. 
When two agents encounter in a trade, their pre-trade wealths 
V, w change into the post-trade wealths v* , w* according to the 
rule 



■ piv + qiw, 



w 



■q2V + p2W. 



(2.3) 



The interaction coefficients pi and qi are non-negative random 
variables. While q\ denotes the fraction of the second agent's 
wealth transferred to the first agent, the difference pi — ^2 is 
the relative gain (or loss) of wealth of the first agent due to 
market risks. We assume that pi and qi have fixed laws, which 
are independent of v and w, and of time. 

In one-dimensional models, the wealth distribution f{t;w) 
of the ensemble coincides with agent density and satisfies 
the associated spatially homogeneous Boltzmann equation of 
Maxwell type ( 11.2b 



d,f + f = Q+{fJ). 



(2.4) 



on the real half line, w > 0. The collisional gain operator 
2+ {t',v), which quantifies the gain of wealth v at time t due to 
binary trades, acts on test functions (f){\v) as 



e+ (/,/)[<?] 



9{w)Q+{fj){w)dw 
((p{v*) + (piw*))fiv)f{w)dvdw, 



(2.5) 



with (•) denoting the expectation with respect to the random 
coefficients pt and qi in ( 12.31 ). The large-time behavior of the 
density is heavily dependent of the evolution of the average 
wealth 



M(f) ::=Mi(f) = / wf(t;w)dw, 



(2.6) 



Conservative models are such that the average wealth of the 
society is conserved with time, M{t) = M, where the value 
of M is finite. In terms of the interaction coefficients, this is 
equivalent to (pi+qi) = {pi +q\) = 1 ■ 

The Boltzmann equation ( 12.4b belongs to the Maxwell type. 
As briefly described in the Introduction, in the Boltzmann 
equation for Maxwell molecules the collision frequency is in- 
dependent of the relative velocity |8], and the loss term in the 
collision operator is linear. This introduces a great simplifi- 
cation, that allows to use most of the well established tech- 
niques developed for the three-dimensional spatially homoge- 
neous Boltzmann equation for Maxwell molecules in the field 
of wealth redistribution. 



A. Pointwise conservative models 

The first explicit description of a binary wealth exchange 
model dates back to Angle |1] (although the intimate rela- 
tion to statistical mechanics was only described about one 
decade later i20ll26ll ): in each binary interaction, winner and 
loser are randomly chosen, and the loser yields a random frac- 
tion of his wealth to the winner. From here, Chakraborti and 
Chakrabarti Ill3ll developed the class of strictly conservative 
exchange models, which preserve the total wealth in each in- 
dividual trade. 



In its most basic version, the microscopic interaction is de- 
termined by one single parameter X G (0,1), which is the 
global saving propensity. In interactions, each agent keeps the 
corresponding fraction of his pre-trade wealth, while the rest 
(1 — A)(v + w) is equally shared among the two trade partners. 



= Av + i(l-A)(v + w), 



= Aw+-(l-A)(v + w). 



w 



v + w. 



(2.7) 



(2.8) 

In result, all agents become equally rich eventually. Non- 
deterministic variants of the model have been proposed, where 
the amount (1 — A)(v + w) is not equally shared, but in a 
stochastic way: 

v* = Av+e(l-A)(v+w), w* = Aw+(l-£)(l-A)(v+w), 

(2.9) 
with a random variable e € (0, 1). 



B. Conservative in the mean models 

Cordier et al. 1I19I1 have introduced the CPT model, which 
breaks with the paradigm of strict conservation. The idea is 
that wealth changes hands for a specific reason: one agent in- 
tends to invest his wealth in some asset, property etc. in pos- 
session of his trade partner Typically, such investments bear 
some risk, and either provide the buyer with some additional 
wealth, or lead to the loss of wealth in a non-deterministic 
way. An easy realization of this idea [29] consists in coupling 
the previously discussed rules ( 12.8b with some risky invest- 
ment that yields an immediate gain or loss proportional to the 
current wealth of the investing agent, 

y* = (A + 77i)v+(l-A)w, w* = (a + 772)w+(1-A)v, 

(2.10) 
The coefficients rji,rj2 are random parameters, which are in- 
dependent of V and w, and distributed so that always v*,w* > 
0, i.e. T7i, 772 > —A. For centered 7],, 

(v*+w*) = (l + (77i))v+(l + (j72))w = v + w, (2.11) 

implying conservation of the average wealth. Various spe- 
cific choices for the 77, have been discussed ||29I1 . The easiest 
one leading to interesting results is Tj,- = ±/i, where each sign 
comes with probability 1/2. The factor jj. G (0,A) should be 
understood as the intrinsic risk of the market: it quantifies the 
fraction of wealth agents are willing to gamble on. 



III. BOLTZMANN EQUILIBRIA 

In conservative markets, where M{t) = M, the details of the 
binary trade determine the profile of the steady wealth distri- 
bution. We introduce the characteristic function 

S(^) = ^(L(/^' +??))-!' (3.12) 



which is convex in s > 0, with S(0) — 1. Also, S(l) =0 be- 
cause of the conservation property ( 12.1 It . The results from 
mili] imply the following. Unless S(s) > for all s > 0, any 
solution f{t;w) tends to a steady wealth distribution Poo{w) = 
/oo(w), which depends on the initial wealth distribution only 
through the conserved mean wealth M > 0. Moreover, exactly 
one of the following is true: 

(PT) if S(a) =0 for some a > 1, then P„{w) has a Pareto 
tail of index a; 

(ST) ifS{s) < for all s > 1, then Poo(w) has a slim tail; 

(DD) if S(a) = for some < a < 1, then P^{w) = 5o(w), a 
Dirac Delta at w = 0. 

In case (PT), exactly the moments Ms{t) with s> a blow up 
as f ^ oo, giving rise to a Pareto tail of index a. We emphasize 
that /(f ; w) possesses finite moments of all orders at any finite 
time. The Pareto tail forms in the limit t ^ °o. 

In case (ST), all moments converge to limits Ms{t) -^ M*, 
so the tail is slim. 

In case (DD), all moments Ms{t) with s > I blow up. The 
underlying process is a separation of wealth as time increases: 
while more and more agents become extremely poor, fewer 
and fewer agents possess essentially the entire wealth of the 
society. 



A. The legacy of Maxwell molecules 

The main tool to obtain the aforementioned results is the 
use of the Fourier transform. This idea, which goes back to the 



seminal work of Bobylev JTlIlt], is well-suited to treat collision 
kernels of Maxwellian type. In particular, the Fourier repre- 
sentation is particularly adapted to the use of various Fourier 
metrics. 



According to the collision rule (I2.3l l. the transformed gain 
term reads 

(3.13) 
The initial conditions turn into 

/o(0)-l and /o'(0) = /M. 

Hence, the Boltzmann equation ( 12.4b can be rewritten as 



dt 



-f{t;^) = \{Kp^^)Kq^^] 



-fiPimqi^))- 

(3.14) 

Equation ( 13.141 1 can be easily treated from a mathematical 
point of view owing to the well-known techniques introduced 
so far to study the Boltzmann equation for Maxwell molecules 
and its caricatures, mainly Kac equation 1271 . 
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